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Abstract. In this paper, we study periodic linear systems on periodic time scales which 
include not only discrete and continuous dynamical systems but also systems with a 
mixture of discrete and continuous parts (e.g. hybrid dynamical systems). We develop 
a comprehensive Floquet theory including Lyapunov transformations and their various 
stability preserving properties, a unified Floquet theorem which establishes a canonical 
Floquet decomposition on time scales in terms of the generalized exponential function, 
and use these results to study homogeneous as well as nonhomogeneous periodic prob- 
lems. Furthermore, we explore the connection between Floquet multipliers and Floquet 
exponents via monodromy operators in this general setting and establish a spectral 
mapping theorem on time scales. Finally, we show this unified Floquet theory has the 
desirable property that stability characteristics of the original system can be determined 
via placement of an associated (but time varying) system's poles in the complex plane. 
We include several examples to show the utility of this theory. 
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1. Introduction 

It is widely known that the stability characteristics of a nonautonomous p-periodic 
linear system of differential or difference equations can be characterized completely by a 
corresponding autonomous linear system of differential or difference equations by a periodic 
Lyapunov transformation of variables [8] I23|, 13 Xj . One application of Lyapunov transfor- 
mations has been in generating different state variable descriptions of linear time invariant 
systems because different state variable descriptions correspond to different and perhaps 
more advantageous points of view in determining the system's output characteristics. This 
is useful in signals and systems applications for the simple fact that different descriptions 
of state variables allow usage of linear algebra to design and study the internal structure of 
a system. Having the ability to change the internal structure without changing the input- 
output behavior of the system is useful for identifying implementations of these systems 
that optimize some performance criteria that may not be directly related to input-output 
behavior, such as numerical effects of round-off error in a computer-based systems imple- 
mentation. For example, using a transformation of variables on a discrete time non-diagonal 
2x2 system, one can obtain a diagonal system matrix which separates the state update 
into two decoupled first-order difference equations, and, because of its simple structure, this 
form of the state variable description is very useful for analyzing the system's properties 

Dm 

Without question, the study of periodic systems in general and Floquet theory in par- 
ticular has been central to the differential equations theorist for some time. Researchers 
have explored these topics for ordinary differential equations HU US] [22] [28] [29] [34j [36] , 
partial differential equations [7J [9] US] [25] , differential-algebraic equations p31 [26] , and dis- 
crete dynamical systems [5J (S3] [35] . Certainly [57] is a landmark paper in the area. Not 
surprisingly, Floquet theory has wide ranging effects, including extensions from time vary- 
ing linear systems to time varying nonlinear systems of differential equations of the form 
x 1 = f(t,x), where f(t,x) is smooth and w-periodic in t. The paper by Shi [33j ensures 
the global existence of solutions and proves that this system is topologically equivalent to 
an autonomous system y' = g(y) via an w-periodic transformation of variables. The the- 
ory has also been extended by Weikard [35] to nonautonomous linear systems of the form 
z = A{x)z where A : C — > C™ x " is an w-periodic function in the complex variable x, whose 
solutions are meromorphic. With the assumption that A(x) is bounded at the ends of the 
period strip, it is shown that there exists a fundamental solution of the form P{x)e Jx with 
a certain constant matrix J and function P which is rational in the variable e 27rM /". 

In a relatively recent paper by Teplinskii and Tcplinskii [35] , Lyapunov transformations 
and discrete Floquet theory are extended to countable systems in ioo(N, M). It is proved that 
the countable time varying system can be represented by a countable time invariant system 
provided its finite-dimensional approximations can also be represented by time invariant 
systems. 

Lyapunov transformations and Floquet theory have also been used to analyze the sta- 
bility characteristics of quasilinear systems with periodically varying parameters. In 1994, 
Pandiyan and Sinha [28j introduced a new technique for the investigation of these systems 
based on the fact that all quasilinear periodic systems can be replaced by similar systems 
whose linear parts are time-invariant, via the well known Lyapunov-Floquet transformation. 

In the paper by Demir [13j , the equivalent of Floquet theory is developed for periodi- 
cally time- varying systems of linear DAEs: 4f{C(t)x) + G(t)x — where the n x n matrices 
C(-) (not full rank in general) and G(-) are periodic. This result is developed for a direct 
application to oscillators which are ubiquitous in physical systems: gravitational, mechan- 
ical, biological, etc., and especially in electronic and optical ones. For example, in radio 
frequency communication systems, they are used for frequency translation of information 
signals and for channel selection. Oscillators are also present in digital electronic systems 
which require a time reference, i.e., a clock signal, in order to synchronize operations. All 
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physical systems, and in particular electronic ones, are corrupted by undesired perturbations 
such as random thermal noise, substrate and supply noise, etc. Hence, signals generated 
by practical oscillators are not perfectly periodic. This performance limiting factor in elec- 
tronic systems is also analyzed in |13) and a theory for nonlinear perturbation analysis of 
oscillators described by a system of DAEs is developed. 

The intent of this paper is to extend the current results of continuous and discrete 
Floquet theory to the more general case of an arbitrary periodic time scale (a special case of 
a measure chain [181 I19[ L defined as any closed subset of K. In particular, the main result 
shows that there exists a not necessarily constant nxn matrix R(t) such that en (to +p, to) — 
&A(to + P,to), where $>A(t,to) is the transition matrix for the p-periodic system x A (t) — 
A(t)x(t), and that the transition matrix can be represented by the product of a p-periodic 
Lyapunov transformation matrix and a time scale matrix exponential, i.e. &A(t, to) = 
L(t)en(t, to), which is known as the Floquet decomposition of the transition matrix $>A(t, to). 
Even though the matrix R(t) is in general time varying, because of its construction, one can 
still analyze the stability of the original p-periodic system by the eigenvalues of R(t), just 
as in the familiar (and special) cases of discrete and continuous Floquet theory. We make a 
note that with this unified Floquet theory, the matrix R(t) does become a constant matrix 
R in the cases that the time scale is R or Z, as one would expect. 

There has been work on generalizing the Floquet decomposition to the time scales 
case in the very nice paper by Ahlbrandt and Ridenhour 3 . However, there are some 
important distinctions between their work and this one. First, Ahlbrandt and Ridenhour 
use a different definition of a periodic time scale. Furthermore and very importantly — 
their Floquet decomposition theorem employs the usual exponential function whereas our 
approach is in terms of the generalized time scale exponential function. Finally, we go on 
to develop a complete Floquet theory including Lyapunov transformations and stability, 
Floquet multipliers, Floquet exponents, and apply this theory to questions of stability of 
periodic linear systems. 

Additionally, this paper develops a solution to the previously open question of the 
existence of a solution matrix R(t) to the matrix equation eu(t, r) = M, where R(t) and 
M are nxn matrices, and M is constant and nonsingular. Whereas the question of a 
generalized time scales logarithm remains open, the solution contained within this paper 
offers a step towards fulfilling this gap in time scales analysis. 

This paper is organized as follows. In Section [2] the generalized Lyapunov transforma- 
tion for time scales is developed and it is shown that the change of variables using the time 
scales version of this transformation preserves the stability properties of the system. Then 
in Section [3J the notion of a periodic time scale is presented and the main theorem, the 
unified and extended version of Floquet's theorem, is introduced for the homogeneous and 
nonhomogeneous cases of a periodic system on a periodic time scale. Three examples are 
given in Section[4]to show that the unified theory of Floquet is functional in the cases T = M, 
T = Z, and more interestingly, when T = Pij. Section [5] introduces the unified theorems for 
Floquet multipliers, Floquet exponents, as well as a generalized spectral mapping theorem 
for time scales. In Section [6l the examples from Section [4] are revisited and the theorems 
introduced in Section [5] are illustrated. We end the paper with the Conclusions, where the 
main ideas are stated and the results of the paper are summarized. In the Appendix, the 
theory of time scales is introduced and necessary definitions and results are stated to keep 
the paper relatively self-contained. An excellent introduction to the subject of time scales 
analysis can be found in Bohner and Peterson's introductory texts [5j [6] . 



2. The Lyapunov Transformation and Stability 

We begin by analyzing the stability preserving property associated with a change of 
variables using a Lyapunov transformation on the regressive time varying linear dynamic 
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system 

x A (t) = A(t)x(t), x(t ) = x . (2.1) 

Definition 2.1. A Lyapunov transformation is an inver tible matrix L(t) G C r 1 d (T, M" xn ) 
with the property that, for some positive 77, p G M, 

||i(t)||<p and detL(t)>r) (2.2) 

for all t G T. 

The two following lemmas can be found in the classic text by A.C. Aitken [J. 

Lemma 2.2. Suppose that A(t) is an n x n matrix such that A~ l {t) exists for all t G T. If 
there exists a constant a > such that ||A _1 (t)|| < a for each t, then there exists a constant 
(3 such that | det A(t)\ > f3 for all t € T. 

Lemma 2.3. Suppose that A(t) is an n x n matrix such that A~ l {t) exists for all t G T. 
Then 

llA m - \det A(t)\ 

for all t € T. 



A consequence of Lemma 12.21 and Lemma 12.31 is that the inverse of a Lyapunov trans- 
formation is also bounded. An equivalent condition to (|2.2p is that there exists a p > such 
that 

\\L(t)\\<p and Hi-^ll <p 

for all t E T. 

Definition 2.4. The time varying linear dynamic equation (|2.1[) is called uniformly stable if 
there exists a finite positive constant 7 such that for any x(to) the corresponding solution 
satisfies 

||z(t)|| <7l|z(fo)ll, t>t . 
Uniform stability can also be characterized using the following theorem. 

Theorem 2.5. The time varying linear dynamic equation (|2.1[) is uniformly stable if and 
only if there exists a 7 > such that the transition matrix $>a satisfies \\<f>A(t,to)\\ < 7 for 
all t > to with t, tg 6 T. 

Definition 2.6. The time varying linear dynamic equation (|2.ip is called uniformly expo- 
nentially stable if there exists finite positive constants 7, A with — A G 1Z + such that for any 
to, x(to) the corresponding solution satisfies 

\\x(t)\\ < \\x(t )\\'ye-x{t,to), t>t . 

Uniform exponential stability can also be characterized using the following theorem. 



Theorem 2.7. The time varying linear dynamic equation (|2.ip is uniformly exponentially 
stable if and only if there exists an A, 7 > with —X G 1Z + such that the transition matrix 
satisfies 

\\3> A {t,to)\\ <7e-x(Mo) 

for all t>U) with t, to G T. 

The last stability definition given uses a uniformity condition to conclude exponential 
stability. 

Definition 2.8. The linear state equation (|2.1[) is defined to be uniformly asymptotically 
stable if it is uniformly stable and given any S > 0, there exists a T > so that for any 
to and x(to), the corresponding solution x(t) satisfies 

\\x(t)\\ < S\\x(t )\\, t>t + T. (2.3) 
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It is noted that the time T that must pass before the norm of the solution satisfies (12.3 
and the constant 8 > is independent of the initial time to- 



Theorem 2.9. The linear state equation (|2.1|) is uniformly exponentially stable if and only 
if it is uniformly asymptotically stable. 



Proof. Suppose that the system (|2.ip is uniformly exponentially stable. This implies that 
there exist constants 7, A > with —A G 1Z + so that ||$^(t,r)|| < 7e_>(i,T) for t > t. 
Clearly, this implies uniform stability. Now, given a 8 > 0, we choose a sufficiently large 
positive constant T > such that to + T S T and e^\{to +T,t ) < -. Then for any to and 
xo, and t > T + to with t, T + to G T, 

\\x(t)\\ = |IM*>*o)*o|| 

< \\$ A (t,t )\\ ||io|| 

< 7e_A(i,io)||a;o|| 
<je-x(to+T,t )\\xo\\ 
<8\\x \\, t>t +T. 

Thus, (|2.ip is uniformly asymptotically stable. 

Now suppose the converse. By definition of uniform asymptotic stability, (12. ip is uni- 
formly stable. Thus, there exists a constant 7 > so that 

\\$A(t,T)\\ < 7 , for all t>r. (2.4) 

Choosing 8 — \, let T be a positive constant so that to +T 1 € T and (|2.3p is satisfied. Given 
a to and letting x be so that ||x a || = 1, we have 

||$ A (i +T,toK|| = \\$ A (t + T,t )\\. 

When xo — x a , the solution x(t) of (|2.ip satisfies 

||z(t +T)|| = \\<S> A (to+T,to)xa\\ = \\$ A (t +T,t )\\ ||x„|| < i||x a ||. 
From this, we obtain 

||*^(*o + r,*o)ll < |- (2-5) 

It can be seen that for any to there exists claimed. Therefore, the above inequal- 

ity (g3]) holds for any t S T. 

Using the bound from the inequalities (|2.4|) and (|2.5p . we have the following set of 
inequalities on intervals in the time scale of the form [r + kT, r + (k + 1)T) T , with arbitrary 
t: 

||M*,i-)II<7, te[r,T + T) T 

\\9 A (t,T)\\ =\\$ A (t,T + T)<Z> A (T + T,T)\\ 

<\\<t> A (t,T + T)\\\\<t> A (T + T,T)\\ 

<\, te[r + T,T + 2T) T 
\\$ A (t,T)\\ = \\t> A (t,T + 2T)<P a (t + 2T,t + T)$ a (t + T,t)\\ 

< ||$ A (t,r+2T)|| \\^ a (t + 2T,t + T)\\ ||$ a (t + T,t)\\ 

<J, te [r + 2T,r + 3T) T . 



In general, for any r g T, we have 

||Mt,r)|| < J, <G [r + fcT,r+(fc+l)T) 
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We now choose the bounds to obtain a decaying exponential bound. We now find a positive 
constant A (with —A € 1Z + ) that satisfies \ < e-\(r + T,r) < e-\(t, r), for t S [r, r + T) T . 
Recall the fact for any positive constant (3 such that —j3 € 1Z + , 



[ lim &(-/3)Ar< /* lim A 7 



t S J to 

Then we must define A such that 

T 



to+T 
to 



ln(-)=/ s hmJ s (-X)Ar. 



It is easy to verify that 

1 1 T 

X 2 



A = lim 

Since T has been established to satisfy (|2.3|) with <5 = ^, for any k £ No and t e [r + fcT, r + (fc + 1)T) T , 
where r + A:T and r + (fc + 1)T are merely the upper bounds on the interval, not necessarily 
elements in the time scale, we now have 

for te[r + kT,T+(k + 1)T) T . 

Then for all t, t e T with t > r, we obtain the correct value of A for the decaying 
exponential bound 

\\$A(t,T)\\ < 7 e-A(i,r). 
Therefore, by Theorem 12.71 we have uniform exponential stability. □ 

Theorem 2.10. Suppose L(t) E C* d (T, M" x "), ™£/i £,(*) invertible for all t e T and A(t) 
is from the dynamic linear system (12. ip . TTien i/ie transition matrix for the system 

Z A (t) = G{t)Z(t), Z{t) = I (2.6) 



where 
is given by 



G(t) = L a ~ (t)A(t)L(t) - L a {t)L A (t) (2.7) 

$ G {t,T) = L- X {t)$ A {t,T)L{T). 

for any t, r € T. 

Proo/. First we see that by definition, G(t) e C rd (T, R nx "). For any r E T, we define 

X(t) =L- 1 (t)$ A (t,r)L(r). (2.8) 

It is obvious that for t = r, X(t) = I. Temporarily rearranging f|2.8[) and differentiating 
L(t)X(t) with respect to i, we obtain [BJ Theorem 5.3 (iv)] 

L A (t)X(t) + L ff (£)X A (t) - $ A (t, r)L(r) - i4(t)* A (t, r)L(r), 

and thus 

L°{t)X A {t) = A{t)^ A {t,r)L{r) - L A (t)X(t) 

= A(t)$ A (t,T)L(T) - L^ifyL-^t^A^Hj) 
= [A(t)-L A (t)L- 1 (t)]$ A (t,T)L( T ). 

Multiplying both sides by L a ~ l (t) and noting ([2~7]) and ([23]) . 

X A (i) = [L a l (t)A(t) - L^ 1 (t)i A (t)L- 1 (t)]$ A (t,r)L(r) 

= G(t)X(t). 
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This is valid for any r € T. Thus, the transition matrix of X A (t) — G(t)X(t) is $ G (i,T) ~ 
L^ 1 (t)&A(t, t)L(t) . Additionally, if the initial value specified in (|2.6|) was not the identity, 
i.e. Z(t ) = Z ^I, then the solution is X(t) = $ G (t, r)Z . □ 

2.1. Preservation of Uniform Stability. 

Theorem 2.11. Suppose that z(t) = L~ 1 (t)x(t) is a Lyapunov transformation. Then the 
system (|2.ip is uniformly stable if and only if 

z A {t) = [L^ 1 (t)A(t)L(t) - (t)L A (t) \ z(t), z(t ) = z (2.9) 

is uniformly stable. 

Proof. Equation (|2.1| and equation (|2.9j) are related by the change of variables z(t) = 
L~ 1 (t)x(t). By Theorem 12. 101 the relationship between the two transition matrices is 

$ G (Mo) = L-\t)<i> A (t,t )L(to). 

Suppose that (12. 1|) is uniformly stable. Then there exists a 7 > such that H^^io)!! < 
7 for all t,to € T with t > to. Then by Lemma l2~3l and Theorem 12. 5[ we have 

||$ G (Mo)|| = WL-^t^A^^L^W 

<\\L-Ht)\\\\*A(t,to)\\\\L(ta)\\ 
lP n 

< — = 7G, 

»7 

for all t,t eT with £ > t . By Theorem |2~51 since ||$ G (Mo)|| < 7g, the system (|2~^|) is 
uniformly stable. The converse is similar. □ 

2.2. Preservation of Uniform Exponential Stability. 

Theorem 2.12. Suppose that z(t) = L _1 (t)a;(t) is a Lyapunov transformation. Then the 
system (|2.1|) is uniformly exponentially stable if and only if 

z A {t)^\L a ' 1 (t)A{t)L(t)~L' J ' 1 (t)L A (t)]z(t), z(t )=z (2.10) 

is uniformly exponentially stable. 

Proof. Equations (|2.1[) and (|2. 10[) are related by the change of variables z(t) = L~ 1 (t)x(t). 
By Theorem I2.10[ the relationship between the two transition matrices is 

$ G (Mo) = L- 1 (i)$A(i,to)i(*o)- 
Suppose that (|2.ip is uniformly exponentially stable. Then there exists an A, 7 > 
with -A € ft+ such that \\§ A {t,t )\\ < 7e_ A (Mo) for all t > t with t,t Q € T. Then by 
Lemma 12.31 and Theorem 12.71 we have 

||$ G (Mo)|| = Hi-^MMo^to)!! 

<||L- 1 (t)||||$ A (*,*o)lll|i(*o)ll 
7p" 

< e_A(t,t ) = 7 G e_A(Mo), 

for all t, to S T with t > to- 

By Theorem 12. 7| since ||$ G (t,< )|| < 7G e -A(Mo), the system (|2.10f) is uniformly expo- 
nentially stable. The converse is similar. □ 

Corollary 2.13. Suppose that z(t) = L- 1 (t)x(t) is a Lyapunov transformation. Then the 
system (|2.ip is uniformly asymptotically stable if and only if 

z A (t) = [l^ 1 (t)A(t)L(t) - L^ 1 {t)L A (t)] z{t), z(t ) = zo 

is uniformly asymptotically stable. 

Proof. The proof follows from Theorem 12.91 □ 
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3. A Unified Floquet Theory 

We want to make a vital distinction that is necessary for comprehension of the notation 
that will be used in the remainder of the paper. When considering solutions to the linear 
dynamic system (|2.ip . just as in the familiar case of T = K, we may write the transition 
matrix differently, depending on properties of the system matrix. 

In the most general case, the solution may always be expressed as x(t) — <&A(t,to)xo, 
where $a is written as 11, Theorem 3.2] 

$ A (t,t )=I + A(t 1 )Ati+ A(n) / v4(t 2 )At 2 At 1 

** io J to «* io 

' A(n) r A(t 2 ) ■ ■ ■ r*" 1 A{ n )An • • • An + • • • , (3.1) 

fo J to J to 

which generalizes the classical derivation using Picard iterates for a first order linear differ- 
ential equation to any time scale. 

When the system matrix commutes with its integral, i.e. 

A{t) J A(t)At = J A(t)At A(t), for all s, t £ [t , oo) T , 

we may write the solution to (|2.ip as x(t) — e A (t,tQ)xo. This type of system is called a 
Lappo-Danilevskii system by the Russian school; see [2]. In this case, the representation of 
the matrix exponential e A is equivalent to the transition matrix $a in (|3-ip . However, the 
property that separates the transition matrix from the matrix exponential (defined as the 
solution to (|2.ip ) is that 

A(t)e A (t,t ) = e A (t,t )A(t), but A(t)$ A (t,t ) + $ A (t,t )A(t). 

Note that if the system matrix is constant, we may again represent the solution with 
the matrix exponential, either using the previous representation or the infinite series [111 
Theorem 4.1] 



oo 



e A (t,t ) = 2_ J A k h k {t,t ). 

Just as on R and /iZ, one of the main properties that distinguishes <&a from ea in the 
context of a solution to a linear dynamic system of the form (12. ip . is that in general $a 
does not commute with A. However, when the solution to (|2.ip may be expressed as e A , it 
is equivalent to the fact that &a and A necessarily commute for all t E T. 

We now state definitions that will be used throughout the paper. 

Definition 3.1. Let p £ [0, oo). Then the time scale T is p-periodic if for all t S T: 

(i) t€T implies t + p e T, 

(ii) n(t)=n(t+p). 

Definition 3.2. A : T -> W ixn is p-periodic if A(t) = A(t + p) for all t G T. 

3.1. The Homogeneous Equation. Let T be a p-periodic time scale. Consider the re- 
gressive time varying linear dynamic initial value problem 

x A (t) = A(t)x(t), x(t ) = x , (3.2) 

where A(t) is p-periodic for all t € T and the time scale T is also p-periodic. Although it 
is not necessary that the period of A(t) and the period of the time scale to be equal, we 
will assume so for simplicity. Furthermore, we will assume that < [i(t) < p and that the 
linear dynamic system and the time scale under discussion are p-periodic, unless references 
to more general systems are made. 

We begin with a definition compulsory for the derivation of the matrix that will be 
introduced in Theorem [ 
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Definition 3.3. Given an n X n nonsingular matrix M and any r G R, wc define the 
principal value of the real power of the matrix M by 



y; T(r + 1) / M-XJ 
^ Q j\T(r-j + l) \ Xi 



where rrii is the multiplicity of the eigenvalue Xi and the Pi are the projection matrices of 
M (See Appendix). 

Remark 3.4. Only the principal values of the matrix functions and the eigenvalues will be 
considered in the remainder of the paper. 

The following theorem is crucial for the development of the unified Floquet decompo- 
sition. 

Theorem 3.5 (Construction of the R Matrix). Given a nonsingular constant matrix M 
and constant p > 0, a solution of the time scale matrix exponential equation e R (to+p, to) = 
M is given by R : T — > C nxn , where 

Mi - I 

R(t) := lim . (3.3) 

Iff has constant graininess on the interval [to, to +p] T , then R(t) is constant. 

Proof. Observe that R(t) // R(t)At = // R(r)ArR(t) for all t,t G T. We can represent 
the transition matrix solution <t>jj(t, to) of the linear dynamic system 

z A (t) =R(t)z(t), z(t ) = z , 

using the matrix exponential. In other words, due to the commutativity of R(t), we have 

•MMo) = e R (t,t ). 

It can be verified by direct calculation that e^(t,to) = R(t)en(t,to) = en(t,to)R(t) for all 
t G T. Using (|3.3p . we obtain 

e R (t,to) =M^. 

To prove this claim, first note that e R (to, to) = M° = I. Delta differentiating, we have 

f+M ( f )-t t-tp 

A/ s ,• M 3 - M~^~ 

e R (t,t ) = lim 

*\n(t) s 

( f(.t) \ 
M~ -I\ , t-t 
lim M p 

s\n(t) s J 

= R{t)e R (t,t ). 



It follows straightforwardly that 

e R (t + p, t ) = M —p— = M. 

□ 

An interesting and useful property of the matrix R constructed in Theorem l3.5l is stated 
in the following corollary. 

Corollary 3.6. The matrices R(t) and M have identical eigenvectors. 
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Proof. Implementing Theorem 19.101 for any of the k <n distinct eigenpairs {Xi,Vi} of M 
we have 

_ s f x _ A 

Mvi = XiVi =>• lim Mpvi = lim A p Vi R(t)vi = lira — - 

i A^ —1 

Thus, the k < n distinct eigenpairs of R(t) are {ji{t),Vi}f =1 , where 7i(t) := lim s \^( t ) — . 

□ 

Lemma 3.7. Suppose that T is a p-periodic time scale and P(t) G 1Z(T, W nxn ) is also 
p-periodic. Then the solution of the dynamic matrix initial value problem 

Z A (t) = P(t)Z(t), Z(t ) = Z , (3.4) 

is unique up to a period p shift. That is, $p(t,to) — &p(t + kp, to + kp), for all t G T and 
k G N . 

Proof. By [6], the unique solution to (|3.4p is Z(t) — <&p(t,to)Zo- Observe 

$$(t,to)Z = P(t)$ P (t,to)Zo and S P (t, t Q )\ t=to Z a = * P (i , t )Z Q = Z . 

Now we show that <J>p(i, to) = $p(t + kp, to + kp), for all £ G T and k G N . We do so by 
observing that $p(< + kp, to + kp)Zo also solves the matrix initial value problem (|3.4|) . We 
see that 

$p t+fcp (t + kp, t + kp)Zo = P(t + kp) $ P (< + kp, t + kp) 

= P(t) $ P (t + kp,t + kp) , and 
$ P (i + kp,t + kp)\ t+ kp=t +kp = $p(* + fcp,to + fep)|t=t 

= $ P (i + kp,t + kp)Z 

= z . 

Since the solution to (|3.4[) is unique, we conclude 

$p(i + fcp, t + fcp) = $ P (<, t ), for all t G T and fc G N . 

□ 

The next theorem is the unified and extended version of Floquet decomposition for 
p-periodic time varying linear dynamic systems. 

Theorem 3.8 (The Unified Floquet Decomposition). The transition matrix for a 
p-periodic A(t) can be written in the form 

® A (t,T) = L(t)e R (t,T)L- 1 (T) forallt,TEY, (3.5) 

where R : T — * C nx " and L(t) G C I 1 d (T, C" xn ) are both p-periodic and invertible at each 
t G T. We refer to (|3.5[) as the Floquet decomposition for $,4. 



Proof. We define the matrix R as in Theorem 13.51 with M := &A(to + P,to)- Using this 
definition, R satisfies the equation 

ep(io + P, h) = ®A(to + P, to). 

Define the matrix L(t) by 

L(t) :^^ A (t,to)e^(t,t Q ). (3.6) 
By definition, L(t) G C£ d (T, C" xn ) and is invertible at each t G T. Now 

$a(Mo) = L(t)e R (t,to), (3.7) 

yields 

*a(«o,*) = e^it^L-^t) = e R (t ,t)L-\t), 
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which, together with (13. 6|) . implies 

$a(*,t) =L(t)e R (t, T )L- 1 {T), 

for all r, t E T. 

We conclude by showing that L(£) is p-periodic. By (|3. 6|) and Lemma [37T1 
L(t + p) = $A(i+P,io)e fl 1 (t +p,to) 

= $ A (t+p, h +p)$A(to + P,to)e R (t ,t +p)e R (t +p,t+p) 
= $> A (t,t )$>A(to +P, t a )e] l 1 (to+p,t )e R (t ,t) 

= L{t). 

□ 

Theorem 3.9. Let $yi(i,£o) = L(t)e R (t,t ) be a Floquet decomposition for <&a- Then 
x(t) = $>A(t,to)x is a solution of the p-periodic system (|3.2p if and only if z(t) = L~ 1 (t)x(t) 
is a solution of the system 

z A (t) = R{t)z(t), z(t Q ) = x . (3.8) 
Proof. Suppose x{t) is a solution to (|3.2[) . Then x(t) = &A(t, t )x = L(t)e R (t, to)xo- Setting 

z(t) := L~ x (t)x(t) = L- 1 (t)L(t)e B (t,t )a;o = e«(t, i )z , 
it follows from the construction of i?(t) that z(i) is a solution of (|3.8p . 

Now suppose z(t) = L~ 1 (t)x(t) is a solution of the system f|3 . 8(1 . Thenz(i) = e_R(t, to)£o- 
Set x(t) := L(t)z{t). Then, 

x(t) = L(t)e R (t,t )x = $ A (t,to)xo, 
so x(t) is a solution of (|3.2p . □ 

Corollary 3.10. XTie solutions of the system Q3.2p are uniformly stable (respectively, uni- 
formly exponentially stable, asymptotically stable) if and only if the solutions of the sys- 
tem (|3.8|) are uniformly stable (respectively, uniformly exponentially stable, asymptotically 
stable). 

Proof. The proof follows from the fact that the systems are related by a Lyapunov change of 
variables and implementation of the appropriate stability preservation theorem in Section [21 

□ 

Theorem 3.11. Given any to G T, there exists an initial state x(to) = Xq 7^ such 
that the solution of (|3 . 2[) is p-periodic if and only if at least one of the eigenvalues of 
e R (t +p,t ) = <& A (to +P,to) is 1. 

Proof. Suppose that given an initial time to with x(to) = xq 0, the solution x(t) is 
p-periodic. By Theorem 13. 8[ there exists a Floquet decomposition of x given by 

x(t) = <f> A (t,to)x = L(t)e R (t,t )L- 1 (t )x . 

Furthermore, 

x(t+p) = L(t + p)e R (t + p,t Q )L^ 1 (t )x = L(t)e R (t + p,t )L^ 1 (t a )x . 
Since x(t) = x(t + p) and L(t) — L(t + p) for each t G T, we have 
e R (t, t )L' 1 (to)x = e R (t + p, t )L' 1 (to)x , 

which implies 

e fl (t,i )i" 1 (io)a;o = e R (t+p,t +p)e R (t Q + p,t Q )L~ 1 (t )x . 
Since e R (t +p,t +p) = e R (t, t ), 

e R (t,t )L~ 1 (t )x Q = e R (t,t )e R (t + p,t )L~ 1 (t )x , 
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and thus 

L~ 1 (t )x a = e R (t +p,to)L~ 1 (to)xo. 

Since L~ 1 (to)xo ^ 0, we see that L~ 1 (to)xo is an eigenvector of the matrix e R (to +p,to) 
corresponding to an eigenvalue of 1. 

Now suppose 1 is an eigenvalue of e R (to +p, to) with corresponding eigenvector zq. Then 
zq is real- valued and nonzero. For any to € T, z(t) = e R (t,to)zo is p-periodic. Since 1 is an 
eigenvalue of e R (to+p, to) with corresponding eigenvector zo and e R (t+p, to+p) = e R (t, to), 

z(t+p) = e R {t +p, t )z 

= e R {t + p, i + p)e R {t a + p, t )z 
= e R (t+p,t +p)z 
= e R (t,t )z 
= z[t). 

Using the Floquet decomposition from Theorem 13.81 and setting xq ■— L(io) z o- we obtain 
the nontrivial solution of (|3 . 2|) . Then 

x(t) = <S> A (t,to)x = L^e^tojL^^xo = L(t)e R (t : t Q )z = L(t)z(t), 

which is p-periodic since L(t) and z(t) are p-periodic. □ 

3.2. The Nonhomogeneous Equation. We now consider the nonhomogeneous uniformly 
regressive time varying linear dynamic initial value problem 

x A (t) =A(t)x(t)+f(t), x(t )=x , (3.9) 

where A e lZ(T,W lxn ), f e C prd (T, M" xl ) n ^(T,R" xl ), and both are p-periodic for all 

tel. 

Lemma 3.12. A solution x(t) of (13. 9p is p-periodic if and only if x(to + p) = x{to). 

Proof. Suppose that x{t) is p-periodic. Then by definition of a periodic function, x{to +p) = 
x(t ). 

Now suppose that there exists a solution of (|3.9p such that x(to + p) — x(to). De- 
fine z(t) = x(t + p) — x(t). By assumption and construction of z(t), we have z(to) = 0. 
Furthermore, 

z A (t) = [A(t+p)x(t+p) + f(t+p)} - [A(t)x(t) + f(t)} 
= A(t) [x(t+p) - x(t)} 
= A(t)z(t). 

By uniqueness of solutions, we see that z(t) = for all t e T. Thus, x(t) = x(t + p) for all 

teT. □ 

The next theorem uses Lemma 13.121 to develop criteria for the existence of p-periodic 
solutions for any p-periodic vector- valued function f(t). 

Theorem 3.13. For all to G T and for all p-periodic f(t), there exists an initial state 
x(to) = xq such that the solution of (|3.9p is p-periodic if and only if there does not exist a 
nonzero z{to) = zo and to £ T such that the homogeneous initial value problem 

z A (t) = A(t)z(t), z(to)=z , (3.10) 

(where A(t) is p-periodic) has a p-periodic solution. 

Proof. The solution of (|3.9p is given by 

x(t) = t )x + f $ A (t, (t(t))/(t)At. 
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From Lemma T3. 121 x(t) is p-periodic if and only if x(to) = x(to + p) which is equivalent to 

r U,+P 

[I-<f> A (t +p,t )]xo 



[° *A(to+P,tr(r))/(T)Ar. (3.11) 
J to 



By Theorem 13.111 we must show that this algebraic equation has a solution for any initial 
condition x(to) and any p-periodic f(t) if and only if e#(to + p, to) has no eigenvalues equal 
to one. 

First, suppose en(t\ + p,t\) = <&A{t\ + P,t\), for some t\ G T, and suppose there are 
no eigenvalues equal to one. This is equivalent to 

det [I-Q A (h +p,h)] ^0. 

Since $a is p-periodic and invertible, we obtain the equivalent statement 

^ det [$ A (to +p,h + p) (I-<P A (h + p,t 1 ))$ A (t 1 ,t )] 

= det [$ A (io +p,h +p)$a(*i,*o) ~ M*o + p,*o)] , (3.12) 

Since $a(^o + f>) = ^(^0) *i) , (|3 . 1 2|) is equivalent to the invertibility of [I — $A(to + 

p, to)], for any to G T. Thus, for any initial time to G T and for any p-periodic f(t), (|3.11[) 
has a solution and that solution is 

/•to+p 

Xo = [I-$ A (to+P,to)]~ 1 / $ A (t +p,(T(T))f(T)AT. 



Now suppose (|3.1ip has a solution for every to and every p-periodic f(t). Given an 
arbitrary to G T, corresponding to any fo, we define a regressive p-periodic vector valued 
function f(t) G C prd (T, R nxl ) by 

f(t) = $ A {a{t),t +p)fo, te[t ,t +p) T , (3.13) 

extending this to the entire time scale T using the periodicity. 
By construction of fit), 

rto+P fto+P 

/ <f> A (to+p,a(T))f(T)AT= / f Ar=pfo. 

J to Jto 

Thus, (|3.1ip becomes 

[I-<f> A (t +p,t )]x =pf . (3.14) 

For any f(t) that is constructed as in (|3.13[) . and thus for any corresponding / , (|3.14[) has 
a solution for xo by assumption. Therefore, 

det [I-3> A (t +p,t )]^0. 

Hence, en(t + p,t ) = ^ A (t + p, to) has no eigenvalue of equal to 1. By Theorem 13. Ill 
(|3 . 10[) has no periodic solution. □ 

4. Examples 

4.1. Discrete Time Example. Consider the time scale T = Z and the regressive (on Z) 
time varying matrix 

_i 2+(-iy 
-l 



A(t) 



2 

2+(-l)' 
2 



which have periods of 1 and 2, respectively. The transition matrix for the homogeneous 
periodic linear system of difference equations 



AX(t) 



, 2+(-l)' 
2 

2+(-l) { _j 



X(t), (4.1) 
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is (from straightforward calculation) given by 



*A(f,0) 



1 



2 t+i 



+ (-V3Y (V3) t+1 + (-73) t+1 



As in Theorem 13.81 consider the equation 



e fl (2,0) = $ A (2,0) = ^3 



(v^) 2 + (-v^) 2 (V3) 3 + (-V3) 3 
(v'S) 3 + (-v'S) 3 (v'S) 2 + (-V3) 2 



which simplifies to 

e«(2,0) = (I + R) 2 = 
Then in the Floquet decomposition is given by 

1 



Bit) := lim — $^(2, 0) = - 1 = 



2 x 





V5 



1 



(4.2) 



1 

2*+i 

1 " 
2 



(v^) 4 + (-V3) 4 (V3) t+1 + (-V3) t+1 

(V3) t+l + (-V3)* +l (Vsy + (-V3Y 





• V3\-t 



(4.3) 



which is constant (as expected) since \x = 1 on Z. Furthermore, the 2-periodic matrix L(t) 
is given by 

L(t) :=$ A (t,0)e- 1 (t,0) 

= $a(*,0) + 

-Vi 

_ l + (-l) 4 v/3+(-l)*(-V3) 
V3+(-l) 4 (-V3) 1 + (-1)* 

Using i? and L(t) above, it is straightforward to verify that indeed (|3.7p . and equiva- 
lently, (|3.5p both hold here. 

4.2. Continuous Time Example. Consider the time scale T = R and the time varying 
matrix 

= \ ~m n 
v ' sin(t) 

which has period 2tt. The transition matrix for the homogeneous periodic linear system of 
differential equations 

-1 
sini 



X(t) 



X(t) 



(4.4) 



1 e *(cosi+sini) 

2 2 



is given by 

M*,o) : 

As in Theorem I3.8[ consider the equation 

e K (2^,0) -$ a (2tt,0) 
Then i2(t) in the Floquet decomposition is given by 

R(t) := Urn -L U^O)^ -/) = J-Log$ A (27r,0) 
which is constant (as expected) since /i e on R. Hence, 



-1 



(4.5) 



en(t, 0) = e 



Hi 



and thus 



-i?t 



- — 1 

2 1 
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Furthermore, the 27r-periodic matrix L(t) is given by 

L(t) = $ A (t,0)e« 1 (*.0) 
= <P A (t,0)e- Rt 

a-* 



1 e ^cosi+sini) 

2 2 

1 ' 

_1 cos t+sin t i 

'"> o J- 



1 



e" u 

2 2 



(4.6) 



Using i? and L(i) above, it is straightforward to verify that indeed (|3.7|) . and equivalently, 
(1331) both hold here. 



4.3. Hybrid Example. This example highlights the scope and generality of Theorem 
on domains with nonconstant graininess. 

Consider the time scale T = Pi.i which has a period of 2 and the regressive (on Pi.i) 
time varying matrix 



A(t) 



g(t) i 

-3 



g(t) = -3 + sin(27rf), 



with period 1. The transition matrix for the homogeneous periodic linear system of dynamic 
equations 

" <?(*) i 

-3 



X(t), 



(4.7) 



is given by 



e fl (i,0) C(t) 
e_ 3 (t,0) 



where (7(t) := / e s (t, a-(r))e_ 3 (r, 0)A? 



Then i?(i) in the Floquet decomposition is given by 

1 

~ W) 



-2e- 3 C(2) 
-2e~ 3 









2 


1 












1 





(4.8) 



which is nonconstant (as expected). Hence 

e R (t,0) -- 



-2e- 3 C(2) 
-2e" 3 



and thus, 



-2e" 3 (7(2) 
-2e" 3 



Furthermore, the 2-periodic matrix L(t) is given by 
L(i)=* A (t,0)ea 1 (t,0) 
" e g {t,0) C(t) 



e_ 3 (t,0) 



2e- 3 (7(2) 
-2e" 3 



(4.9) 



which is obviously p-periodic on P14. Thus, the Floquet decomposition of the transition 
matrix is ^(i, 0) = L(t)eR(t,0). 
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5. Floquet Multipliers and Floquet Exponents 

Suppose that $ A (t, to) is the transition matrix and <fr(t) is the fundamental matrix at 
t = t (i.e. $(r) = /) for the system (13.21) . Then we can write any fundamental matrix &(t) 
as 

*(t) = *(i)tf(r) or *(t) = * A (t,to)*(*o). 

Definition 5.1. Let £0 G M n be a nonzero vector and ^f(t) be any fundamental matrix for 
the system (|3.2|t . The vector solution of the system with initial condition a; (to) = ^0 is given 
by $ A (t,t )x . The operator A/ : M™ -> R" given by 

M(x Q ) := <f> A {t Q +p,t Q )x = ^(to+p^-^^xo, 

is called a monodromy operator. The eigenvalues of the monodromy operator are called the 
Floquet (or characteristic) multipliers of the system (|3.2j) . 



The following theorem establishes that characteristic multipliers are nonzero complex 
numbers intrinsic to the periodic system — they do not depend on the choice of the funda- 
mental matrix. This result is analogous to the theorem dealing with the eigenvalues and 
invertibility of monodromy operators in [8] , which can also be referenced for proof. 

Theorem 5.2. The following statements hold for the system (|3.2[) . 

(1) Every monodromy operator is invertible. In particular, every characteristic multi- 
plier is nonzero. 

(2) // Mi and M2 are monodromy operators, then they have the same eigenvalues. In 
particular, there are exactly n characteristic multipliers, counting multiplicities. 

With the Floquet normal form $a(Mo) = *i(f)*^ x (to) = L(*)e fl (Mo)£~ X (*o) of the 
transition matrix for the system (|3.2[) on one hand, and the monodromy operator represen- 
tation 

M(x ) = <S> A (t + p,t )x = #i(t +p)^i 1 (t )x , 
on the other, together we conclude 

<5> A {t Q +P,to) = 9i{to +p)*r 1 (M = L(t )e R (t + p,^)^ 1 (t Q ). 

Thus, the Floquet (or characteristic) multipliers of the system are the eigenvalues of the ma- 
trix enito+p, to). The (possibly complex) scalar function j(t) is a Floquet (or characteristic) 
exponent of the p-periodic system (|3.2I) if A is a Floquet multiplier and e 7 (to +p, to) = A. 

Theorem 15.31 will help us answer the question of whether or not the eigenvalues of 
the matrix R(t) in the Floquet decomposition & A (t,to) = L(t)eu(t,to) are in fact Floquet 
exponents. 

Theorem 5.3 (Spectral Mapping Theorem for Time Scales). Suppose that R(t) is 
an n x n matrix as in Theorem \3.5[ with eigenvalues 71 (t), ...,7 n (i), repeated according to 
multiplicities. Then Ji(t), 7„(t) are the eigenvalues of R k (t) and the eigenvalues of ex 

Proof. By induction for the dimension n, we start by stating that the theorem is valid for 
lxl matrices. Suppose that it is true for all (n — 1) x (n — 1) matrices. For each fixed 
t G T, take 71 (t) and let v ^ denote a corresponding eigenvector such that R(t)v — ji(t)v. 
Let ei, e„ denote the usual basis of C n . There exists a nonsingular matrix S such that 
Sv = ei. Thus we have SR(t)S~ 1 ei = 7i(t)ei, and the matrix 5i?(t)5 _1 has the block 
form 

SR ^ S ^ = [ 7 (<) R*( t ) 

The matrix SR k (t)S~ 1 has the same block form, only with block diagonal elements 7^ (t) and 
R k (t). Clearly, the eigenvalues of this block matrix are 7* (t) together with the eigenvalues 
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of R k (t). By induction, the eigenvalues of R k (t) are the fcth powers of the eigenvalues of 
R(t). This proves the first statement of the theorem. 

Since we chose the matrix S so that SR^S' 1 is block diagonal, by construction, the 
matrix eg^g-i is also block diagonal. We see that esRS- 1 has block diagonal form, with 
block diagonal elements e 7l and e^. Using induction, it follows that the eigenvalues of 
are e 72 , e 7n . Thus, the eigenvalues of egng-i — Se^S -1 are e 7l , e ln . □ 



We now know that the eigenvalues of the matrix en(to+p, to) arc the Floquet multipliers 
and the eigenvalues of R(t) are Floquet exponents. However, in Theorem 15.51 we will see 
that although the Floquet exponents are the eigenvalues of the matrix R(t), they are not 
unique. We first introduce the definition of a purely imaginary number on an arbitrary time 
scale. 

Definition 5.4. Let — £ < w < £. The Hilger purely imaginary number ilo is defined by 
iu) = - — t=- -. For z G Cft, we have that zlm^(z) G Also, when h = 0, iui = iuj. 

Theorem 5.5 (Nonuniqueness of Floquet Exponents). Suppose j(t) G 1Z is a (possibly 
complex) Floquet exponent, A is the corresponding characteristic multiplier of the p-periodic 
system (|3.2|) such that e 1 (to + P, to) = A, and T is a p-periodic time scale. Then r y(t) 
is also a Floquet exponent for all fcgZ. 



Proof. Observe that for any k G Z and £q 6 T, 



e 7 e?2ni= (*o +P,to) = e~f(to +p, to)eo 2 „ k (t +p,t ) 



e 1 (t +p,t )exp 

e 7 (to +p,to) exp 

e 7 (t +p,< )exp 
e 7 (i +p,t )exp 
e 7 (i + p,i )exp 



to+P Log(l + fj,(r)t 



In 



M(r) 



27rfc \ 

^Ai 



in 



At 



in 



io+PLog(l + MT r * ) 
ta+p Log(e i27r,CAl(T)/p ) 

, W) 

to+p i2irkfi{T)/p 

ta+p i2irk . 
At 



An 



P 



e 7 (t +p,t )e 
e 7 (io +P,*o)< 



i2irk 



□ 



The next lemma will illustrate the periodic nature of a special exponential function 
which will be used in proving the nonuniqueness of Floquet exponents in Theorem 15.71 

Lemma 5.6. Let T be a p-periodic time scale and k G Z. Then the functions eo 27rk (t,to) 
and e (f , to) are p-periodic. 
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Proof. Let t e T. Then 

e°2^c +P,to) = exp 



i2nk(t + p — to) 



( i2nk(t — to) \ / i2nkp 
— exp exp 1 

VP/ 

/i2nk(t - t ) 
= exp 

V p 

Therefore, eo 27rk (t,to) is a p-periodic function. The fact that e o 27rk (t,to) is p-periodic 
follows easily from the identity [5] 

1 

P 

□ 

We now show that for any Floquet exponent, there exists a Floquet decomposition such 
that the Floquet exponent if an eigenvalue of the associated matrix R(t). 

Theorem 5.7 (Special Floquet Decompositions). Ifj(t) is a Floquet exponent for the 
system (13. 2|) and <f>A(t,to) is the associated transition matrix, then there exists a Floquet 
decomposition of the form ^^(i, to) = L(t)eu(t,to) where "fit) is an eigenvalue of R(t). 

Proof. Consider the Floquet decomposition ^(i, to) = L(t)e R (t,to). By definition of the 
characteristic exponents, there is a characteristic multiplier A such that A = e 7 (p + to, to), 
and, by Theorem 15.31 there is an eigenvalue j(t) of R(t) such that e 7 (p + to, to) — A. Also, 
by Theorem 15.51 there is some integer k such that j(t) = j(t) ® ?2^. 
Set 



. . ~ , . o27rfc 

R(t) := R(t) e i I, 

P 

L{t) :=£(t)e S2ai (t,io). 

p 

By this definition it is implied that Rit) = R(t) i^^-I. Then j(t) is an eigenvalue of R(t), 
L{t) is a p-periodic function, and 

L(t)e R it,t Q ) = L(t)eo 3 „ h Jt, t )e R (t,t ) = Ljt)e a ^ k ImR jt,to) = L{t)e R (t,t ). 

p p 

It follows that $A(t,to) = Lit)enit,to) is another Floquet decomposition where j(t) is an 
eigenvalue of Rit). □ 

The following theorem classifies the types of solutions that can arise with periodic 
systems. 

Theorem 5.8. If X is a characteristic multiplier of the p-periodic system (13. 2[) and e 7 (fo + 
p, to) — A for some to G T. then there exists a ipossibly complex) nontrivial solution of the 
form 

x{t) = e 7 it,t )qit) 

where q is a p-periodic function. Moreover, for this solution x(i + p) = Ax(t). 

Proof. Let $>A(t,to) be the transition matrix for (|3.2|) . By Theorem 15. 7\ there is a Floquet 
decomposition <&yt(i, to) = L(t)eR(t, to) such that j(t) is an eigenvalue of R(t). There exists 
a vector v ^ such that R(t)v = 7(t)v. It follows that en(t, t )v — e T (t, t )v, and therefore 
the solution x{t) := &A(t, to) v can be represented in the form 

x(t) = Lit)e R (t,t )v = &y(t,to)L(t)v. 
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The solution required by the first part of the theorem is obtained by defining q(t) := L(t)v. 
The second part of the theorem follows from 

x(t + p) = e 7 (f + p, t )q(t + p) 

= e 7 (f + p, t + p)e 1 (t + p, t )q(t) 
= e 7 (f + p,f )e 7 (f + p, t + p)q(t) 
= e 7 (f + p 1 t )e 1 (t 1 t )L(t)v 
= e 7 (f + p,t )x(t) 
= Xx(t). 

□ 

The next result is motivated by [El Theorem 2.53] and is a direct consequence of The- 
orem l5~ 



Corollary 5.9. Suppose that Ai, .. . , A„ are the Floquet multipliers for the p- periodic system 
(& 

(1) If all the Floquet multipliers have modulus less than one, then the system (|3.2j) is 
exponentially stable. 

(2) If all of the Floquet multipliers have modulus less than or equal to one, then the 
system (|3.2p is stable. 

(3) If at least one of the Floquet multipliers have modulus greater than one, then the 
system (|3.2p is unstable. 

Theorem 15.81 showed that if 7(f) is a Floquet exponent of ()3.2|) . then we can construct 
a nontrivial solution of the form x(t) = e~(t,ta)q(t), where q(t) is p-periodic. In the next 
theorem it is shown that if two characteristic multipliers Ai and A2 of the system (|3.2[) are 
distinct, then as in Theorem l5.8[ we can construct linearly independent solutions x% and xi 

of Era. 



Theorem 5.10. Suppose that Ai, A2 are characteristic multipliers of the p-periodic system 
()3.2j) and 71 (£), 72(f) are Floquet exponents such that e 7i (fo+p, fo) = Aj, i = 1,2. If Ai 7^ A2, 
then there exist p-periodic functions 91(f), 92(f) swc/i f/iaf 

xi(f) = e 7l (f,f )(?i(f) and x 2 (f) = e 72 (f, f )g 2 (f) 
are linearly independent solutions of Q3.2p . 



Proof. As in Theorem 15.71 let $^(f, fo) = L(t)en(t,to) be such that 71(f) is an eigenvalue 
of i?(f) with corresponding (nonzero) eigenvector «i. Since A2 is an eigenvalue of the mon- 
odromy matrix <&A(fo +P, to), by Theorem 1 5 . 31 there is an eigenvalue 7(f) of R(t) such that 
e 7 (fo +P,to) = A 2 = e 72 (f +p,f ). Hence 72(f) = 7(f) © for some k E Z. Also, 
7(f) 7^ 71(f) since Ai 7^ A2. Thus, if i>2 is a nonzero eigenvector of R(t) corresponding to the 
eigenvalue 7(f), then the eigenvectors v\ and vi are linearly independent. 
As in the proof of Theorem 15.81 there are solutions of (|3.2[) of the form 

x\(t) = e 7l (f, t )L(t)vi, x 2 (t) = e 7 (f,f )L(f)u 2 . 

Because xi(fo) = V% and a^2(fo) = 1^2, these solutions are linearly independent. Finally, xi 
can be written as 

%2 (f) = ^ e 7ffi °2 ILil (f,fo)^ [e e ^(t,t )L(t)v2 \ =e 72 (f,f ) (f , t )L(t)v 2 

where q 2 (f) := (t, t )L(t)v 2 . □ 
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6. Examples Revisited 

We now revisit the examples from Section 01 We show for each of the three examples in 
Section 2] that the original R(t) matrices and the corresponding Floquet exponents are not 
unique. We conclude each example with a stability analysis using the theorems developed 
from Section [3l 



6.1. Discrete Time Example. In Section [4.11 the Floquet exponent we found for the 



system was ^ — 1. We show that 7 := — — 1 is also a Floquet exponent, but it is not 



an eigenvalue of the original matrix R in 

Set R:= RQ ilixlkjp = R Q iwl, with k = 1 and p = 2, as in Theorem ETT1 Then 



R(t) 



4-1 

VI _ 

2 













ITT 







e 







1 







ITT 



(fi)e» 



(±f 1 >K 



.VI _1 

2 ± 



Hence, 



e A (t,0) = (I + R)* 





_ VI _ 1 
2 1 



Vj 

2 





2 



and 



(-1)* 

(-iy 



eojt,0) = (-iy 

Next, using L(t) from (|4.3p . set L(i) := L(t)e<> Jt, 0) to obtain 



L(t) = 



_ l + (-l) 4 _ V3 + (-l)'(-V3) 
V3+(-l)*(-V3) l + (-l)' 

1-1)* \ _ l + (-l) 4 _ V3+(-l)*(-V3) 
V3 + (-!)*(- v^3) ! + (-!)* 



(-1)' 
(-l) f 



Thus, 

and so 
L(t)e B (t,0) 



£(i)e A (i, 0) = L{t)e %I {t, 0)e A (t, 0) = L^e^^t, 0) = L(t)e K (t, 0) 



_ 1+ (-1)* _ >/3 + (-l)*(-V3) 
V3 + (-l)*(-V3) ! + (-!)* 

_ l+(-l)* _ >/3 + (-l)*(-V3) 
>/3+(-l)*(-V3) ! + (-!)' 



(-1)* 
(-1)* 



(-*)' 



_ (-1)* [ _!+(-!)' _ V3+(-l)'(-V3) 
2 [ ^+(-l) t (- v / 3) 1 + (-!)* 

= L(tWt,0). 



2 



VT 
2 , 
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Therefore, <&a(£, 0) = L(t)e^(t,0) is another (distinct) Floquet decomposition of the tran- 
sition matrix for A. Moreover, 



7 



:= (V3/2 - 1) e !tt = -V5/2 - 1, 



is a Floquet exponent as well as an eigenvalue of R which corresponds to the Floquet 
multiplier A = |; that is, 

W-i)e^( 2 '°)= e -V3/ 2 -i(2 ; 0) = 3/4. 

In light of Corollary 15.91 solutions of (|4.ip are exponentially stable since the Floquet 
multipliers Ai,2 = 3/4 have modulus less than 1. 



6.2. Continuous Time Example. From Section T4.21 consider R from (14. 5[) . Set R :- 

R i2irlk/p = R — iirl, with k = 1 and p = 27r,as in Theorem 15.71 Then, 



R = 



"-10" 




" i " 




" 


. 1 o_ 




i 




. 1 



so that 



and 



e«(t,0) = e R * 



$a(27t,0) = 



e (-l-i)t 

„-<« «,(-!-<)« 



e" 27r 

1 

2 2 



Next, using L(t) from flUD, set L(i) := L(t)eo (t,0) = L(t)e ut to obtain 

e 4 * 



L(i) = 



1 

_1 _ cos t+sin £ i 



e il 




e*' _ e' t (cost+sint) it 
2 2 e 



Thus, 



and so 



L(t)eg(t,0) = L(t)e m e Rt = L(i) e ( iJ+i ^ = L(i)e ffl , 



1 

]_ cos £+sin £ i 

2 2 

1 

1 cos £+sin t i 

2 2 i 



e 

1 e 



2 2 

e lt 




1 



e 



3 (-i-t)t 

„(-!-<)« 



e %t _ e'^cost+sint) a 
2 2 e 



,(-!-*)* 



= L(t)e 



Therefore, ^(t, 0) = L(i)e^(i,0) is another (distinct) Floquet decomposition of the tran- 
sition matrix for A. Moreover, 71 := — 1 — i, 72 := — i are Floquet exponents as well as 
eigenvalues of R which correspond to the Floquet multipliers Ai = e~ 2lr , A2 = 1, respectively. 
That is, e - 2 *- 2m = er 2 * and e~ 2 ™ = 1. 

In light of Corollarv l5.91 solutions of (I4.4[) are uniformly stable since the Floquet mul- 
tipliers satisfy |A» I < 1, i = 1,2. 
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6.3. Hybrid Example. In Section l4~3l consider R from (|4.8j) which was given by 



R(t) 



-2e~ 3 




C(2) 
-2e~ 3 









2 


" 1 












1 





where — 2e 3 = e_3(2,0) on T. The eigenvalues of R(t) are 

7l (i):= lim s- 1 ((e_ 3 (2 ) 0))-'-l) and j 2 {t) := lim s" 1 ((e_ 3+sin(27rt) (2, 0))" 



!)■ 



Set i? := i? e i2nlk/p = RQ ml, with = 1 and p = 2, as in Theorem 1571 Then, 



and thus, 



-2e" 3 (7(2) 
-2e" 3 











o 




2 


" 1 


" 




ITT 










1 


• 







ITT 



2e~ 3 


(7(2) 


2 





-2e~ 3 




2e~ 3 


(7(2) 


t 
2 





-2e- 3 




2e~ 3 


(7(2) 


t 
2 





-2e" 3 







— iirt 



Next, using L(t) from (|4.9|) . set Z(i) := L(t)e a (t, 0) to obtain 
L(t) 



e 9 (*,0) C(t) 
e_ a (t,0) 



-2e- 3 (7(2) 
-2e" 3 



2 















Thus, 



L(t)e s (t, 0) = L(t)e* nI (t, 0)e & (t, 0) = L(t)e k (t, 0)e^(t, 0) = 0) = L(t)e fl (i, 0), 

and so 



L(t)e B (t,0) = 



e S (*>0) r e s (t, ( r(T)) e _3(r,0)AT 
e_ 3 (t,0) 



-2e- 3 / 2 e 9 (2, ( r(r)) e _3(r,0)AT 
-2e~ 3 



-2e- 3 / 2 es (2, t T(T)) e _ 3 (r,0)Ar 



-2e" 



" e a (t,0) 


r o t e ff (i,a(r)) e _ 3 (T,0)AT 




" -2e" 





e-s(*,0) 







-2e- 3 


/ 2 e 9 (2, ( r(r)) e _3(r,0)Ar " 


t 
2 







-2e~ 3 







" e g (t,0) 


r o t eg (t,a(r)) e _ 3 (r,0)Ar 




" -2e" 





e_ 3 (*,0) 







-2e- 3 


/ 2 e 9 (2, ( r(r)) e _3(r,0)Ar " 


t 
2 







-2e" 3 








-2e 



-3 



-ITTt 

r2 





— ITTt 





' e i7rt 


■ 


e Mt 



-2e" 



= L(i)e 5 (<,0). 

Therefore, ^(i, 0) = L(t)e^(t,0) is another (distinct) Floquet decomposition of the tran- 
sition matrix for A. Moreover, 7(f) := —3 G iir is another Floquet exponent as well as an 
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eigenvalue of R(t) which corresponds to the Floquet multipliers X1.2 



-2e 



-3 



That is, 



senr 



(2,0) = e_ 3 (2,0) = -2e- 



In light of Corollary 15. 9( solutions of (|4.7|1 are exponentially stable since the Floquet 
multipliers satisfy |A»| < 1, i = 1,2. 



7. Floquet Theory, Stability, and Dynamic Eigenpairs 
Consider the p-periodic regressive system 

x A (t) = A{t)x(t), x(t ) = x . 
Recall that the R matrix in the Floquet decomposition of $^ is given by 

R(t) := lim - (§ A (to + P,to)' ~ r 

s\fJ,{t) S \ 



(7.1) 
(7.2) 

(7.3) 



and consider the uniformly regressive system associated with (|7.ip . 

z A (t) = R(t)z(t), z(t ) = x a . 

From Theorem l3.9[ solutions to these two problems are related via z(t) = L _1 (i)x(i) where 
L{t) is the Lyapunov transformation from (|3 . 6[) . 

Given a constant n x n matrix M, let C be the nonsingular matrix that transforms M 
into its Jordan canonical form, 

J := C-HlC = diag [J mi (A x ), . . . , J mfe (A fc )] , 

where k < n, 5Zi=i m « = n i ^ are the eigenvalues of M (where some of them may be equal), 
and J m (A) isamxra Jordan block, 



J m (A) = 



A 1 
A 



1 

A 



(7.4) 



We introduce a definition from [30) which gives a bound on the eigenvalues of a system. Then 
we state a lemma which proves that the system (|7.3j) associated with (|7.ip via a Floquet 
decomposition is necessarily uniformly regressive. 

Definition 7.1. [30] The scalar function 7 : T — > C is uniformly regressive if there exists a 
constant 5 > such that < 5 1 < |1 + |, for all t G T K . 

Lemma 7.2. Each eigenvalue of the matrix R(t) in (|7.3[) is uniformly regressive. 

Proof. We show by direct substitution that the eigenvalues of R(t) are uniformly regressive. 

_ a^— 1 

Using Corollary I3.6[ let JiW '■— ^ m -s\p.(t) ~ l ~ — be any of the k < n distinct eigenvalues of 

R(t). Recall that in this p-periodic setting < /z(t) < p. Suppose that |Aj| > 1. Then, 

xf - 1 



|l + /i(t)7i(*)l = lim 



1 



On the other hand, suppose < |Aj| < 1. Then, 

xf - 1 



|l + p(i)7<(*)l= lim 



1 



= lim |A, P I > 1. 



lim \Xf I > |A, 



Thus <5 1 := min{l, |Ai|, . . . , |Afe|} will suffice as the bound for uniform regressivity. 
The following definition will be used to obtain the main result of this section. 



□ 
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Definition 7.3. Let C M := {z e C : z ^ Given an element t E T K with /i(i) > 0, 

we define the HUger circle as 

H t :={zeC, : |l + */i(t)| < 1}. 
Thus, we define the closed Hilger circle as 

H t :={zeC M :|l + z^(<)l<l}, 
Likewise, when fj,(t) = 0, we define the Hilger circle as 

Hi := {z e C : Re 2 < 0} 

and the closed Hilger circle as 

Hi := {z e C : Rez < 0}. 

Motivated by the work of Wu [37] on R, in [12] we introduced the following. 

Definition 7.4. A nonzero, delta differentiable vector w(t) is said to be a dynamic eigen- 
vector of M(t) associated with the dynamic eigenvalue if the pair satisfies the dynamic 
eigenvalue problem 

w A {t)=M(t)w(t)-£(t)w a {t), teT K . (7.5) 

We call {£(£), w(t)} a dynamic eigenpair. 
The nonzero, delta differentiable vector 

roi(i) := e^(t,to)wi(t), (7.6) 

is the mode vector of M(t) associated with the dynamic eigenpair Wi(t)}. 

The following lemma proves that given any regressive linear dynamic system (not neces- 
sarily periodic), there always exists a set of n dynamic eigenpairs, with linearly independent 
dynamic eigenvectors. 

Lemma 7.5. Given the n x n regressive matrix M, there always exists a set of n dynamic 
eigenpairs with linearly independent dynamic eigenvectors. Each of the eigenpairs satisfies 
the vector dynamic eigenvalue problem (|7.5p associated with M. Furthermore, when the 
n vectors form the columns of W(t), then W(t) satisfies the equivalent matrix dynamic 
eigenvalue problem 

W A (t) = M(t)W(t) - W a (t)E(t), where S(t) := diag[£i(i), . . .,£»(«)]. (7.7) 

Proof. Let {£i(t)}2=i be a set of (not necessarily distinct) regressive functions. Then the 
n x n nonsingular matrix W{t) defined by 

W(t) :=<f> M (t,t Q )e e3 (t,t ), 

has as its columns the associated n linearly independent dynamic eigenvectors {wi(t)}™ =1 . 
By direct substitution into (|7.7[) . the proof is complete. □ 

The next theorem is from [12] . We show that the stability of a linear dynamic system 
can be completely determined by the mode vectors m,, 1 < i < n, as constructed in (|7.6p . 

Theorem 7.6 (Stability Via Modal Vectors). Solutions to the uniformly regressive (but 
not necessarily periodic) time varying linear dynamic system (|2.ip are: 

(i) stable if and only if there exists a 7 > such that every mode vector rrij(t) of A(t) 
satisfies ||mj(i)|| < 7 < 00, t > t , for all 1 < i < n, 

(ii) asymptotically stable if and only if in addition to (i), ||mj(t)|| — > 0, t > to, for all 
1 < i < n, 

(iii) exponentially stable if and only if there exists 7, A > with —A G 7?. + (T,]R) such 
that \\mi(t)\\ < jex(t,to), t > to, for all 1 < i < n. 
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Proof. Let {£,i(t) , Wi(t)}f =1 be a set of n dynamic eigenpairs with linearly independent dy- 
namic eigenvectors associated with the system matrix in (|2.1|L The transition matrix can 
be represented by 

$ A (t, h) = W(t)e B {t, ^W^ito) (7.8) 

where W{t) := [ Wl {t), W 2 (t),--- ,w n (t)] and H(t) := diagfci(i), • • • , £»(*)]■ 
Denoting the matrix W (to) as 

" «i"(*o) 

T 
"2 



W- L (t ) 



V T(t ) 



«»(*o) 



we now have the reciprocal basis of each iWj(to) given by the row vectors vf (to). 
Because S(i) is a diagonal matrix, (|7.8[) can be rewritten as 



M*>*o) = £e C4 (t J f )W r (i)F i W- 1 (i ) 1 



(7.9) 



where Fi :— Sij is n x n. Observing that vf (t)wj(t) — §ij for all t G T, we may rewrite 
as 

^ = VK- 1 (i)[0,...,0, tB j (t),0,... J 0]. (7.10) 
Substituting (|7TTU|) into (JTSJ) yields 



i=l 



The proof concludes easily from (|7.11|) . 



(7.11) 



□ 



We now set the stage for the main result of the section. In the next theorem, we 
will show that given the system matrix R(t) from (|7.3[) . we can choose a set of n linearly 
independent dynamic eigenvectors that have a growth rate that is bounded by a finite sum 
of generalized polynomials. 

Theorem 7.7. Given the set of traditional eigenvalues {"ti(t)}f =1 from the matrix R(t) in 
(|7.2p . let {wi(t)}^ =1 denote the corresponding linearly independent dynamic eigenvectors as 
defined by Lemma \ 7. 5\ Then {-fi(t) , Wi(t)}f =1 is a set of dynamic eigenpairs of R(t) with 
the property that each Wi(t) is bounded by at most a finite sum of constant multiples of 
generalized polynomials. In other words, there exists positive constants Di > such that 

IK(*)|| < A j2 h k {t,t ), (7.12) 

k=0 

where rrii is the dimension of the Jordan block which contains the ith eigenvalue, for all 
1 < % < n. 

Proof. The fact that {■~fi(t),Wi(t)}™ =1 is a set of dynamic eigenpairs of R(t) follows imme- 
diately from Lemma 17.51 We now show each dynamic eigenvector Wi satisfies the bound 
in (|7.12p . For an appropriately chosen nonsingular n x n constant matrix C, we define the 
Jordan form of the constant matrix <&A{to + P, to) by 

J :=C- l § A {to+p,to)C (7.13) 

Jmi (^l) 

JmA^d) 
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where d < n, 5Z i=1 rrii — n, Xi are the eigenvalues of &A(to +P,to) (where some of them 
may be equal), and J m (A) is a m x m Jordan block as defined in (|7.4|) . 
Using the matrix C from (|7.13p . define 



K(t) — C-^R^C 
= C~ 1 lim 

1 



(7.14) 



s V 



= lim -(c-^Aito+p^'C-l) 



7 (*"')■ 



= lim 

where the penultimate equality is via Theorem 19.61 Note that K (t) has the block diagonal 
form 

- K^t) 

K 2 (t) 



K(t) 



where each Ki (t) from (|7. 15|) has the form 



K d (t) 



(7.15) 



Ki{t) 



>'(*) 

A, p -1 



A 



A. p -1 



2!p 



A i " -1 



(J £ffi_ 1) ... (J £Mi_ ( „_ 2))Ai P 



(n-l)!p 
fi-(„-3))A, i 



-(n-2) 



(n-2)!; 



l!p 
tt£*3 
A, p -1 
M (t) 



By hypothesis, the dynamic eigenvalues of i?(i) have been defined as its traditional 

A^ — 1 

eigenvalues 7,(i) = lim s ^ M ( t ) — ^ — , for 1 < i < n. Note that since R(t) and are similar, 
they have the same traditional eigenvalues, with corresponding multiplicities. Furthermore, 
choosing the n dynamic eigenvalues of K(t) to be its n traditional eigenvalues, we claim the 
corresponding dynamic eigenvectors {ui(t)}f =1 are defined to be Ui(t) ;= C~ 1 w i (t). 

To prove this, since {"fi(t), Wi(t)}f =1 is a set of dynamic eigenpairs for R(t), by definition 
we have 

wf(t) = R(t)wi(t) - -Yi(t)w?{t), 

for all 1 < i < n. We can now show that {ji(t), Ui(t)}" =1 is a set of dynamic eigenpairs for 
K(t): 

u?(t) = C- 1 w?t(t) (7.16) 
= C- 1 R{t)w i (t)-C- 1 li {t)w° i {t) 
= K^C^Wiit) ~ 7i (i)C- 1 <(t) 
= K{t)u i (t)- li {t)u1{i), 

for all 1 < i < n. 

We now show that each of the Wj(t) is bounded by a finite sum of constant multiples of 
generalized polynomials. 

Since each of the n dynamic eigenpairs {7i(t), satisfy the dynamic eigenvalue 

problem (|7.16p . we can now find the structure of each of the dynamic eigenvectors. Choose 
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the ith block of K(t), with dimension mm x rm. We must then solve the m s x m, linear 
dynamic system 



v*[t)=Ki{t)v{t) 



n J_ M«)-p) 

U Xi P 2X' z p 2 



1 



(M(t)-p)-(M(*)-p("-2)) 
(n-1)! A" _1 p n_1 

(mW-p)-(p-'W-p("-3)) 

(n-2)! A^ 2 p»- 2 



1 

AiP 





Kt), (7-17) 



where -Kj(i) := Ki(t) G 7i(i)l. Since there are m« linearly independent solutions to (|7.17p . 
we will denote each of the rrij x 1 solutions by Vij(t), where i corresponds to the zth block 



matrix Ki(t) of K(t) and j = 1, . . . , mi. For 1 < i < d, define U := X)s=o ms ' w ith mo 
The form for an arbitrary u x 1 column vector u^+j, with 1 < j < mi, is 

«£+,■(*) = [0, . . • , 0, . . . , 0, 0, . . . , 0, . . . , Olixn- 



0. 



(7.18) 



miH hm ; _i 



mt+iH hmj 



By combining all n vector solutions from (|7.16[) , the solution to the equivalent nx n matrix 
dynamic equation 

U A (t) = K{t)U{t) - U a (t)T(t), 
where T(t) := diag [71 (t), ■ ■ ■ , 7n(*)L will have the form 



U(t) 



Ul ,...,U mi ,..., U (E *-i mfc ), ■ ■ • , U( Ei=i mfc ), ■ ■ • , U( E * =i mfc )_!, U„ 
Wl.l «1,2 ■ • ■ Ul.mi 



«1,1 



«l,mi -1 



Vd,2 ■ ■ ■ Vd,m d 
Vd,l '• Vd,m d -1 

Vd,l 



m d Xm d J nX n 



For brevity, we focus on the m, elements contained in our m, x 1 solution vector j/ with 
the understanding that it will be embedded into an n x 1 vector of the form (|7.18p . in the 
appropriate place, so that it only acts on the corresponding block matrix Ki(t) of K{t) 
in (|7.16p . There are rrii linearly independent solutions of (|7.17p and they have the form 

:= [«i,m«(*),0, • . -,0]^. xl 

^i,2(*) := [«i,m«-l(*),Ui, m< (t),0,...,0]^. xl 



[vi,2(t),v i)3 (t), . . . ,Vi tmi -i(t),v iimi (t),0] miXl 

[Vi,l(t),Vi 7 2(t),Vi 73 (t), . . . 1 Vi irn ^i(t),V t ^ mi (t)]^ n , xl , 
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where the scalar functions Vij are constructed by back solving based on the form of Ki in 
(|7.17|) . The associated dynamic equations are 



vtm,(t)=0 (7-19) 



\iP 



1 



A / \ — P ( \ 1 ( \ 

ZAjP AiP 
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-p)0(*) 
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3)P) 


(/*(*) 
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1 

AiP 



/iff) — p . . 1 

Ui,m 4 -l(*) H 1- -^TTl~ v ^\ t ) + T - u »,3(*j 

2AjP Aip 



a /.\ \y\" y \y\" ■ ■ ■ \y\°J ^-w /.\ 
"WW = 7- 77r7^7=r— i Vi,mM 



Vi, mi -i(t) h h 2 «i,3(t) + - — Vi^yt). 

l\P AiP 



The solutions to (f7TT9)) are 

U»,m, (*) = 1 



Vi, mi -l(t) 



f l 1 

/ T V iymi (T)AT 

J t MP 



f u(r) — p . /" 1 . . . 

Vi, mi -2(t) = / — ■r r5 -t)j, mi t)At + / - — D^i-i t At 
y* 2Afp2 7 tQ A 4 p 

/■* (/i(r)-p)(/i(r)-2p) f' nW-P f , A 



3!A 3 p3 J tQ 2A?p2 

/■* 1 

/ t — Ui, mi _2(T)Ar 



4 (/x(r) - p)(/x(r) - 2p) . . . (p(r) - (to, - 3)p) 

Ut, m< -i(r)Ar 



w <.a(*)=/ ; ,,,.^-2 m - 2 u iimj (r)Ar 

t (mj - 2)!A i p m > ^ 

* (p(r) - p)(p(r) - 2p) . . . Oi(r) - (m, - 4)p) 
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(to 2 -3)!A™^V 1 *- 3 



J t (m 4 - l)!A i p m « 1 

* (Kt) - 2p) . . . (/i(r) - (to, - 3)p) 

Ui, mi _i(TjAr 

to 



(to, - 2)!A? li ~V H - 2 
V(r)-P / ,\ a , f* 1 



+ --- + L^^ v ^ mT+ L-p v ^ T - 
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We can explicitly bound each Vij which yields an explicit bound on each u^+j, for all 
1 < j <rrii and for all £i, with 1 < i < d. Recall /j,(t) < fi max < P for all tel. There exists 
constants Bij, i = 1, . . . , d, and j = 1, . . . , m i; such that 

|w mi (i)| = 1 < Bi^„ u h (t,t ) = Bi^ 1l%1 
\v m ,-i(t)\ < [ ^-v mi (T)AT = M^iM < B itmt -ihi(t,to), 



\v mi -2(t)\ < 



H{r)-p 



2A?p 2 



v mi (r) 



At 



to 



■T—V mi -i(T) 

Kp 



At 



< 



A?P 2 



\v mi -a(t)\ < [ 

<-f 

Jto 



hi(t,t ) h 2 (t,t ) 
~ 2\\p + \\p> 

2 

< B iiTni _2 hj(t, tp), 
J'=l 

1 ([i(t)-p)(p(t)-2 P ) 



At 



3!A 3 p 3 



v mi (t) 



At + 
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Jtn 
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v mi -i(r) 



At + 
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At + 
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3!Afp 3 
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M < B i>2 ^2 hjtt^o), 

— 1 

M < B iA ^2 hj(t,t ). 

For each 1 < i < d, set /?, := max. (= i i ... iTOi {.B, ; j}. Then, for 1 < i < d and j = 1,.. ., m i} we 
have 



\\u h+j (t)\\ <fo hk(t,t ). 



fc=0 



.'SO 



DACUNHA AND DAVIS 



Finally, recall that Wi — Cu L using C from (j7. 14|) and define Di :— ||C||/3j, for all 
1 < i < n. Then we have the bound 

m;-l rrii — 1 

\\w i {t)\\ = \\Cu i {t)\\<\\C\\p i h k (t,to)=Di h k (t,to), 

k=0 k=0 

which proves the second claim in the theorem. □ 

To complete the setup for the main theorem, we state a necessary result from [3"TJ which 
characterizes the growth rates of generalized polynomials. 

Lemma 7.8. Given a regressive scalar function A : T — > C, \£ C^ d (T, C) on a time scale 
with bounded graininess and unbounded above with the property 

3T e T : < - inf te[Ti00)T [Re M (A(t))] , 

it holds that 

lim h k (t,t )e x (t,T) = 0, /orreT, fceNo. 

t — >oo 

We are now in position to prove the main stability result. The result will show that 
given a p-periodic time varying linear dynamic system (|7.1[) , the traditional eigenvalues of 
the associated time varying linear dynamic system (|7.3|) (via the Floquet decomposition of 
$a) completely determine the stability characteristics of the original system. 

Theorem 7.9 (Floquet Stability Theorem). Given the p-periodic system (|7. 3|) with 
eigenvalues {"fi(t)}^ =1 , we have the following properties of the solutions to the original p- 
periodic system (|7.ip : 

(i) if Re^i{i) < for all i = 1, . . . ,n, then all the system (|7.1[) is exponentially stable; 

(ii) if i?e M 7i(t) = for all i — l,...,n, then in the case when only simple elementary 
divisors correspond to "fi{t), the system (|7.1|) is stable. The system (17. 1|) is unstable, 
growing at rates of generalized polynomials of t if among the elementary divisors 
there are multiple ones; 

(hi) if Re^i(t) > for some i = 1, . . . ,n, then the system (17. ip is unstable. 

Proof. Let <&A(t, to) be the transition matrix for the system (|7.1| . Define the matrix R(t) as 
in (|7.2|) . and note that by Theorem 13.81 (|7.3|) is the associated system through the Floquet 
decomposition of Given the traditional eigenvalues {7i(i)}™ = i of R{t), we define a 

set of dynamic eigenpairs of R(t) by {ji{t), u>i(i)}™ =1 as in Lemma [731 By Theorem 17.71 
the dynamic eigenvectors Wi(t) are bounded by a finite sum of generalized polynomials as 
in (|7.12|) . Employing Theorem 17.61 we can write the transition matrix of (|7.3|) as 

n 

ZR.{t, t ) = ^ m »( t )^ T ( t o), 
i=l 

for all 1 < i < n, where the vectors mi(t) and vj {to) are defined as in Theorem 17.61 
Case (i): Applying Lemma [7T51 for each 1 < i < n, we have 

rm — 1 

lim ||mi(i)|| < D l V h k (t, t )|e 7i (t, t )\ < C £ e Rc [li]S)e (t,t Q ) < C e e-\(t,t ) = 0, 

k=0 

where e > is chosen so that 7i(t)©e G H t , for all t E T, C e := max t > to Di Y^k=o 1 ^fc(^ *o)ee e (t, £o), 
and A > with Re M [7i(t)] ® e < -A < 0, for alH = 1, . . . , n and tel. By Theorem 1761 
(hi), the system (|7.3| is exponentially stable. 

The solution to (|7.1|) and the solution to (|7.3jl are related by a Lyapunov transformation, 
namely L(t) :— &A(t, io)e^. 1 (t, ^o)- By Theorem 13.91 and Corollary 13 . 10j the solution to 
the system (|7.1I) is exponentially stable if and only if the solution to the system (|7.3|) is 
exponentially stable. 
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Case (ii): Suppose Re^ [71(f)] = for all i = l,...,n with only simple elementary 
divisors corresponding to each Ji(t). Then each Jordan block corresponding to each 7i(f) is 
lxl, which implies that the mode vectors have the form 

rm(t) = Ae 7i (f,f ). 

Thus, 

lim ||mj(f)|| < lim ^e Re u, i (t, t ) = Pi < 00. 

By Theorem 17.61 (ii), the system (|7.3j) is stable. Thus, since (|7.ip is related to ()7.3|1 by a 
Lyapunov transformation, (|7.1jl is stable. 

Now suppose that for some i the elementary divisor corresponding to 7$ (f ) has multiplic- 
ity n > 1. Then, the elements in the corresponding mode vector is made up of generalized 
polynomials. Thus, since there is no decaying exponential function to control the growth 
of these polynomials, solutions grow at the rate of the polynomials. Therefore, (|7.3|) is 
unstable, implying by the Lyapunov transformation that (|7.1|) is also unstable. 

Case (iii): Suppose that for some 1 < i < n, we have Re^ [71(f)] > 0. This implies that 

lim ||e fl (f,fo)|| = 00, 

t — >oo 

thus implying by the Lyapunov transformation 

lim ||$ A (f,f )|| = 00. 

t— >oo 

□ 



The consequences of Theorem 17.91 should be underscored at this point: from the con- 
struction of the corresponding system matrix R(t), we can deduce the stability of (|7.3[) 
solely on the placement of the traditional eigenvalues in the complex plane. In general, this 
is not true of time varying systems. However, we have shown that for this unified Floquet 
theory, the associated (but in general time varying) matrix R(t) yields stability character- 
istics through simple pole placement in the complex plane — exactly like constant systems 
and Jordan reducible time varying systems on K or hZ. Moreover, the matrix R(t) is in fact 
constant if and only if the domain of the underlying dynamical system has constant graini- 
ness, e.g. R or hZ. The upshot of course is that this Floquet theory extends to any domain 
that is a p-periodic closed subset of R, including those that have nonconstant discrete points 
as well as mixed continuous and discrete intervals. 

A key observation about this unified Floquet theory is that the associated system 
matrix R(t) yields the stability characteristics of the original system by placement of the 
traditional eigenvalues alone — although it accomplishes this via an associated time varying 
systems rather than with a time invariant system. 

Remark 7.10. We can think of an eigenvalue 7(f) = lim^s^t) ~* of R(t) as constant 
with respect to its "placement" in the instantaneous stability region corresponding to the el- 
ement f in the p-periodic time scale T. In other words, the exponential bound is independent 
of /i(f): 

e 7 (t, fo) = exp 

"/"*,. t ( A* - A 
= exp / urn Log 1 + s 

[Jt »\,Mt) y s j 



lim 



Log(l + sj(t)) At 
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8. Conclusions 

A unified and extended Floquet theory has been developed, which includes a canonical 
Floquet decomposition in terms of a generalized matrix exponential function. The notion 
of time varying Floquet exponents and their relationship to constant Floquet multipliers 
has been introduced. We presented a closed-form answer to the open problem of finding 
a solution matrix R(t) to the equation en(t,T) = M on an arbitrary time scale T, where 
R{t) and M are n x n matrices, and M is constant and nonsingular. The development of 
this Floquet theory on both homogeneous and nonhomogeneous hybrid dynamical systems 
was completed by developing Lyapunov transformations and a stability analysis of the orig- 
inal system in terms of a corresponding (but not necessarily autonomous) system that, by 
its construction, yields stability characteristics by traditional eigenvalue placement in the 
complex plane. 

9. Appendix: Real Powers of a Matrix 

We introduce the technical preliminaries required to define the real power of a matrix, 
which is of paramount importance in the derivation of the unified Floquet theory. The 
following can be found in [20l [21] . 

Definition 9.1. Given an n x n invertible matrix M with elementary divisors {(A — 
Ai) mi }f =1 , we let the characteristic polynomial be the function p(X). We define the polyno- 
mial Oj(A) implicitly via the expansion 

1 _ ^ a»(A) 

The polynomial bi(X) is defined by omitting the factor (A — Xi) mi from the characteristic 
polynomial p(X), that is 

bi(X) : 1I' :A A ^ • 

Then the ith projection matrix is defined as 

Pi(X) := ai (X)bi(X). (9.1) 
A property of the set of the projection matrices is the following. By definition, 

1 = E f A -A A U p(A) =5X A ) 6 iW =E p i( A )- ( 9 - 2 ) 

3 = 1 K 3 ' 3=1 j=l 

When the scalar value A is replaced by the n x n nonsingular matrix M, (|9.2p becomes 

k 

I = J2Pj(M)- (9-3) 

We now state and prove three propositions concerning the projection matrices. 

Proposition 9.2. Given an n x n nonsingular matrix M as in Definition Iff. 11 we have 
Pi(M)(M- Xil) mi = 0. 

Proof. Observe 

k 

Pi{M){M-XiI) mi = a l (M)6,(Af)(M-A l /) m - = a t (M) JJ(M-A J -I) m * = a t (M)p(M) = 0. 

□ 



Proposition 9.3. The set {P^M)}^ is orthogonal. That is, P i (M)P j (M) = SijPi{M). 
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Proof. Suppose i ^ j. Then 

Pi(M)Pj(M) = Pi (M)a,j (M)bj (M) 

= Pi(M) aj (M) J[(M - X k I) m - 

= Pi(M) aj (M){M - Xil) mi J] (M - A fe /r fc 

= Pi(M)(M - \J) m >a,(M) J] (M - A fe /)" lfc 
= 0. 

To conclude the proof, using the property (|9.3| and multiplying on each side by the ith 
projection matrix, we have 

k 

P t {M) =Y J Pi{M)P j {M) = Pi(M)Pi(M). 

3 = 1 

□ 

Proposition 9.4. Given an n x n nonsingular matrix M as in Definition \9.1\ with corre- 
sponding Jordan canonical form J = CAIC^ 1 , for some nonsingular matrix C , we have 

CP^NDC- 1 = Pi(J), 

for all 1 < i < k. 

Proof. For any 1 < i < k, by definition of the projection matrix (|9.ip . we see that Pi(M) 
is the product of two polynomials of the matrix M. Further, by definition of ai and bi, it 
follows that Ca l (M)C' 1 = a t {J) and Cb^M)^ 1 = b t (J). Thus, 

CP l (M)C- 1 = Ca,i(M)bi{M)C- x = C a^M)^ 1 Cb^M)^ 1 = ai {J)bi{J) = P l (J). 

□ 



We are now in the position to state the definition of the principal value of the real 
power of a nonsingular n x n matrix. 



Definition 9.5. Given an n x n nonsingular matrix M as in Definition 19. II and any r € 

we define the real power of the matrix M by 



— 1 

E 

3=0 



T(r + 1) fM-Xjy 



j\T(r-j + l)\ Xi 



The next theorem gives a relationship between the real power of a matrix and the real 
power of the corresponding Jordan canonical form. 

Theorem 9.6. Given annxn nonsingular matrix M as in Definition \9. li vith corresponding 
Jordan canonical form J = CMC^ 1 , for an appropriately defined nonsingular matrix C, 
and any r € R, we have 

CM r C- x = J T . 
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Proof. The result follows from a direct application of Proposition 19.41 and Definition 19.51 
Observe, 



CM r C- 1 =cY^Pi(M)\l 



mi — 1 

E 

3=0 



r(r + l) (M-XilV 



j!r(r-j + l)V A 



c- 1 



= ^CP l {M)C- 1 C\ r l 



k 

= Y J Pi{J)K 

i=l 

= J r . 



E; 



T(r + 1) (M-Xil 



j\r(r-j + l)\ X t 



771 ^ — 1 

E 

J=0 



T(r + 1) 



j!r(r-j + l)V A 



J — AjJ 



□ 



Remark 9.7. Theorem 19.61 is a generalization from the familiar matrix property that 
CM k C- 1 = J k , for all k e N . 

Remark 9.8. Given any nonsingular nxn matrix M as in Dcfjnition l9.il we have Pi(M)(M— 
Xil) mi = 0, for i = 1, . . . , k and mi + • • • + mjt = n. However, since Pi(M)(M — Aj7) is 
nilpotent by Proposition 19.21 and Proposition [931 there exists an integer rii < mi such that 
Pi(M)(M - XJ)"*- 1 ^ 0, but Pi{M){M-XiI) ni = 0. Then P,(M)(M - A,;/)"* is a minimal 
polynomial of M. 

Theorem 9.9. j4 nonzero column vector v r of Pi(M)(M — X;J) ni ~ r is a generalized eigen- 
vector of rank r, for 1 < r < n,, associated with the eigenvalue A;. A regular eigenvector 
has rank r = 1 . 

Proof. For r = 1, we have 

fc 

MPi(M)(M - A,/)"'- 1 = £ P i (M) [A,/ + (M - Xjl)] P t (M)(M - XJ)^- 1 

i=i 

= A i P i (M)(M-A^) n '- 1 , 

where we use the result from Proposition 19.31 which states that Pi(M)Pj(M) = 6ijPi(M), 
for all 1 < i,j < k. Thus, since Pi(M)(M — Xil) 71 ^ 1 has rank one, we choose a nonzero 
ii x 1 column vector vi from this matrix and obtain 

Mvi = A,ui. 

When r > 1, note that 

P(M)(M - Ai/) ni ~ r = (M - XiI)Pi[M){M - A^) n '- r - 1 ' 

which implies 

iV = (Af - A i /)w r+ i. 

□ 



The following theorem will show that for any invertible nxn matrix M, with d eigenpairs 
{Xi,Vi}f =1 , and any r <E R, the matrix M r has the d eigenpairs {A^,i>i}^ =1 . 



Theorem 9.10. Given an invertible nxn matrix M as in Definition \9.1\ with d eigenpairs 
{Xi,Vi}f =1 , and any r > 0, for each eigenpair {Xi,Vi} of M, we have M r vi — X^Vi. 
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Proof. Let 1 < s < d and X s be an eigenvalue of the matrix M. Then, noting that (M 
A s J) ns is a minimal polynomial of M, using Proposition 19.21 and Remark 19.81 we have 

M r P s {M){M~ XJY'- 1 

nii — 1 



r(r + l) (M-Xil 



E j! T(r-j 



P S (M)[M - X s I) n ° 



p -mk E j!r r ( ^7li) A ^ (M - A - /)J+ " s 



= A s r F s (M)(M-A s /)"- 1 , 

for each nonzero column vector v s of P S (M)(M-A S /)" S_1 . The matrix P s (M)(M-X s I) ne ~ 1 
has rank 1 and we can take any nonzero column vector v s from this matrix as an eigenvector 
corresponding the eigenvalue X s of the matrix M. This vector is a linear combination of all 
of the eigenvectors associated with the eigenvalue A s . Thus, M r v s — X r s v s . □ 
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